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Abstract - Using Brownian Dynamics (BD) simulations we investigate the non-equilibrium struc¬ 
ture formation of a two-dimensional (2D) binary system of dipolar colloids propelling in opposite 
directions. Despite of a pronounced tendency for chain formation, the system displays a transition 
towards a laned state reminiscent of lane formation in systems with isotropic repulsive interac¬ 
tions. However, the anisotropic dipolar interactions induce novel features: First, the lanes have 
themselves a complex internal structure characterized by chains or clusters. Second, laning occurs 
only in a window of interaction strengths. We interprete our findings by a phase separation process 
and simple force balance arguments. 


Introduction. — Lane formation is a protoype of a 
non-equilibrium self-organization process, where an orig¬ 
inally homogenous mixture of particles (or other types 
of ’’agents”) moving in opposite directions segregates 
into macroscopic lanes composed of different species. 
This ubiquitous phenomenon occurs, e.g., in driven bi¬ 
nary mixtures of colloidal particles nHa and migrating 
macroions [J, in binary plasmas 0E3, but also in ’’self- 
propelling” systems with aligned velocities such as bacte¬ 
ria in channels [7] and humans in pedestrian zones [5]. In 
particular, studies of charged colloids have revealed many 
fundamental aspects of laning such as the impact of den¬ 
sity 0, the role of hydrodyamics pD], the accompanying 
microscopic dynamics (particularly, the so-called dynami¬ 
cal locking) [2 , and the impact of anisotropic friction El- 

All of these models involve isotropic interactions be¬ 
tween the colloidal particles. However, recent progress 
in colloidal chemistry has generated a variety of complex, 
anisotropic particles (see, e.g., 1T21fRl ). which can perform 
controlled translational motion in an external field EES], 
like simple charged colloids, but display complex self- 
assembly behaviour already in equilibrium OS El- In 
driven ensembles of such particles one may therefore ex¬ 
pect a wealth of new phenomena induced by the interplay 
between self-assembly, on the one hand, and dynamical 
self-organization processes such as laning, on the other 
hand. The consequences are so far, only poorly under¬ 
stood, contrary to the widely studied case of active (self- 


propelled) particles ( [T8H20] ). 

In the present letter we discuss, based on particle-based 
Brownian Dynamics (BD) simulations, a prototype of 
field-propelled complex colloids, that is, spheres with dipo¬ 
lar interactions. Our model is inspired by real systems of 
metallodielectric ’’Janus” spheres with two dielectric parts 
acquiring different (induced) dipole moments in an exter¬ 
nal electric field [21] . Previous experimental studies in a 
quasi-two dimensional (2D) setup with an in-plane field 
have shown that these particles perform straight motion 
perpendicular to the field El: with two possible directions 
depending on the orientation of the hemispheres. Ensem¬ 
bles of such particles therefore should exhibit lane forma¬ 
tion, with the additional feature of strongly anisotropic 
interactions favoring the formation of staggered chains. 

Based on a suitable model system, we find indeed sev¬ 
eral new phenomena: First, despite a pronounced ten¬ 
dency to aggregate into chains perpendicular to the driv¬ 
ing force, the dipolar interactions can induce laning at den¬ 
sities where purely repulsive systems are mixed. At larger 
interaction strengths, however, we observe a breakdown of 
lane formation. Moreover, by comparing the lane forma¬ 
tion with that in a simpler system governed by isotropic, 
attractive Lennard-Jones (LJ) interactions, we show that 
the main laning mechanism is, in fact, the effective (angle- 
averaged) attraction between the dipoles. Indeed, for both 
systems the onset of laning roughly occurs at coupling 
strengths related to equilibrium phase separation. Finally, 


P-1 





F. Kogler et al. 



Fig. 1: (Color online) Potential energy ’landscape’ due to 
the total dipolar interaction (see color code) for two DDPs 
of diameter a (solid circles) with (a) Si = Sj = — 1 and (b) 
Si = 1, Sj = —1 in a top-view perspective for (p 1 )* = 1.58. 
The particle dipole moments are indicated by arrows with a 
white [black] head for /A 1 ' [/x <2, [. The white circle around the 
centered particle indicates the excluded area. 


lanes disappear when the thermal energy exceeds the at¬ 
traction (corrected by the drive) between two particles. 

Theoretical models and methods. — The exper¬ 
imental particles motivating this work are gold-patched 
dielectric spheres confined between two glass plates, form¬ 
ing a quasi-2D geometry [15] . Application of an in-plane 
AC electric field aligns the plane between gold patched and 
dielectric part along the field and induces dipole moments 
in both, the gold and the dielectric part of the particles, 
with the gold’s dipole being significantly larger due to its 
larger polarisability. At low frequencies these two induced 
dipoles have the same orientation (along the field). Placed 
into an aqueous solution, one observes m a spontaneous 
motion of each particle in the direction orthogonal to the 
field and away from the particle’s gold patch, the under¬ 
lying mechanism being asymmetric flow of solvent charges 
(induced-charge electrophoresis) [15l|22] ■ 

To examine the collective behaviour of these driven par¬ 
ticles, we perform BD simulations in a 2D quadratic cell 
of size L 2 with periodic boundary conditions. The cell 
contains N = 800 spherical particles with equal diameter 
a defining a unit length. Particle positions are denoted 
by Vi = Xie x + yie v (i = 1,iV) with the unit vectors 
e x and e y defining a 2D coordinate frame. The exter¬ 
nal, electric field E ext points along the y-axis. Sterical 
interactions are modeled by a soft-sphere (SS) potential 
Uss{ r ij) = 4ess (( a / r ij) 12 — (v/rij) 6 + 1/4) truncated 
at rF = 2 s a, where r,j = \rj — 771 is the particle distance. 
The strength of repulsion is set to t* ss = tss/ksT = 10 
where ks is Boltzmann’s constant and T is the temper¬ 
ature. To mimic the impact of the propulsion force we 
randomly assign to each particle a fixed vector Sj = Sie x 
which points either to the ’right’ (s = 1) or to the ’left’ 
(s = —1) with probability 0.5, respectively. Thereby a 
random fifty-fifty mixture of two different particle ’species’ 
is created. Each particle is then subject to a constant but 
orientation-dependent driving force f di = fd^i- Changes 
of particle orientation (sj) are not considered here. 

To incorporate the dipolar interactions we assume that 


each particle i bears two point dipole moments /r-b 
( 2 ) 

and ; whose orientation is fixed along the external 
field, i.e. = /i| Q ^e y , a = 1,2. Moreover, the 

two dipole-moments are shifted out of the particle cen¬ 
ter by 5^ = 5si and 6^ = ~5si with S = 0.25er. 
The values of are hold constant (that is, we neglect 
variations of the local field) and set differently in order 
to mimic the strong asymmetry of the metallo-dielectric 
janus-particles. Specifically, we choose fj,^ = 2/Lt- 1 ^ and 
set the dimensionless dipole strength to /j* = ji) = 


/yJa^ksT. Test simulations revealed that the sys¬ 
tem’s behaviour is only weakly sensitive to the ratio 
/U 2 /Mi, whereas changing the shift parameter 6 has a 
strong effect. Here we stick to a value of S already used 
in Ref. [T7] for the case of oppositely oriented dipoles. 

(a) 

The distance between two arbitrary dipole moments ' 
and /i$ 7) °f different particles is given by r®'* = Tj + 

SJ. — (rq + Sg.). The dipolar interaction between the 
two of them is calculated via the 3D point dipole po- 

*i, Sj ) =^“ ) -MS 7) '- a ' Y ^ 3 


tential U dip (r iy , n[ a) , ^ 


, 0-7 

SiS. 


r 3 - 


3 (/A (q) • (Mj 7) • r s7sj)( r VBj) 5 > and the total dipo¬ 

lar interaction between two particles is given by the sum 

UDD^ij j Si, Sj) = Xya,7=l Udip(f'ij , , Sj ). We 

use the standard 2D Ewald summation method to han¬ 
dle the long-range character of the dipole interactions |23 : . 
Combining Udd and the formerly introduced SS potential 
defines the ’double-dipole particle’ (DDP) model Uddp = 
Udd + Uss- The angle dependency for distances ry > cr 
is illustrated in fig. [T] It is seen that two identical DDPs 
(si = Sj) prefer a head-to-tail configuration similar to 
’’single-dipole” particles with one dipole in their center. 
In contrast, for DDPs with s,; ^ Sj, the most attractive 
configurations are staggered in the sense that one DDP is 
shifted towards the gold-patched ’back’ of the other parti¬ 
cle. The corresponding minimum energy Up^(si ^ Sj) is 
roughly 4 kgT larger than in the head-to-tail condiguration 
Ujjd^($ i = Sj). Finally, for reasons discussed later we also 
consider a driven LJ fluid characterized by the pair po¬ 
tential Ulj = 4e((cr/rij ) 12 — (cr/rij ) 6 with cut-off at 2.5 ct 
and dimensionless interaction strength e* = e/ksT. For 
all model systems the overdamped BD equations of motion 
are given by yrq = VU x (ij)+f s di i+Ci where 7 is the 

friction constant, Ux{ij) is a specific pair potential (with 
X = SS, DDP, or LJ), and is a Gaussian noise vector 
which acts on particle i and fulfills the relation {£,) = 0 
and (Ci(t)Cj(t 0) = 2y kBT5ijS(t — t'). Hydrodynamic in¬ 
teractions (HI) are neglected, in accordance with earlier 
studies [10: revealing that HI do not alter lane formation 
qualitatively. The BD equations are solved via the Euler 
scheme [54] with an integration stepwidth At = 10 _ 5 rb, 
where Tb is a brownian timescale defined by tj, = tr^y/fcaT. 


Target quantities. — To quantify the degree of lan- 
ing we employ the laning order parameter $ defined in [3;. 
Perfectly laned states correspond to $ = 1. However, the 
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Fig. 2: (Color online) Laning order parameter as function of 
driving force for (a) the DDP system ( p* = 1.58) and (b) the 
SS system (etj s = 10 ksT) at different densities. 


system can visually appear as being laned already at much 
smaller values of $ (e.g., $ ~ 0.2), especially at high den¬ 
sities. To characterize the local structure in 2-direction we 
calculate the radial distribution function between particles 
of the same type 

N N 

9 S ( X ) = l a: *il)®( cr — |j/*jl)®( s i s j)) (!) 

^ i=l 

with Xij (yij ) being the x (y)-component of r^, and 8 and 
0 are the delta- and the Heavyside step function, respec¬ 
tively. By normalization, g s {x) decays to 2 for x —> oo in 
a single-species fluid and to 1 in a completly mixed binary 
fluid. By interchanging x and y in eq. © we additionally 
define the correlation function g s (y) in y-direction. 

Results. — In fig. [Ua) we plot the laning order pa¬ 
rameter of DDP systems for different densities p* = 
pa 2 as functions of the dimensionless driving force f* { = 
f d a/k B T. At low driving forces (f d < 10) the order 
parameter is essentially zero for all densities considered, 
reflecting the absence of lanes. The corresponding lo¬ 
cal structure at /^ = 0 and p* = 0.2 is illustrated by 
the simulation snapshot shown in fig- Eta). One observes 
the formation of ’staggered’ chains, consistent with the 
most attractive pair configurations illustrated in fig.[]Jb), 
and in qualitative agreement with experiments of met- 
allodielectric particles under field conditions where self¬ 
propulsion is absent [ 21] , Increasing the driving force 
from zero < 10), these staggered chains are pushed 
against each other and form one large aggregate consist¬ 
ing of both species. Upon further increase of f* : , all of the 
systems considered in fig. 0(a) display laned ’states’ char¬ 
acterised by non-negligible values of <f>. For small densities 
(p* < 0.4) this lane formation occurs gradually and signif¬ 
icant values of $ are reached only at high driving forces 
f d > 50. In contrast, dense systems display a rather steep 
laning ’’transition” at a driving force « 15. A visu¬ 
alisation of exemplary laned states at = 80 and two 
densities is given in fig. 0(b) and (c). 



Fig. 3: (Color online) Snapshots of the DDP system at (a) 
fd = 0, p* = 0.2, (b) f* d = 80, p* = 0.2 ($ « 0.2) and (c) 
f d = 80, p* = 0.5 (4> « 0.75). Dark (bright) particles: s = 1 
(-1), p* = 1.58. 


To elucidate the impact of dipolar forces on the func¬ 
tions 4>(/J) we plot in fig- 0(b) corresponding data for the 
driven SS system. At low densities p* < 0.4 the systems 
behave similar to their DDP counterparts. However, at 
high densities (p* > 0.4) there is no lane formation in 
the (purely repulsive) SS system, in striking contrast to 
the behaviour of the DDP systems. We will come back 
to this point below. Here we first consider the internal 
structure of the lanes in the driven DDP systems. Visual 
inspection of figs. [31(b) and (c) suggests that lanes consist 
of clusters of particles of the same type. At low densi¬ 
ties (see fig. [31(a)) we find mostly short chains oriented 
in y-direction, which reminds of the equilibrium structure 
of single dipole particles in external fields [25]. Clearly, 
such chains are absent in simpler (isotropic) lane-forming 
systems. At higher densities (see fig. 131(b)) the particles 
tend to form clusters along the driving force. Their inter¬ 
nal structure somewhat resembles that in a (hexagonal) 
crystal; in fact, the corresponding translational order pa¬ 
rameter [26] is larger than that observed at low f d . For 
additional quantitative characterization we consider the 
correlation functions g s (y) and g s (x) plotted in fig. [4] For 
the lowest density considered (p* = 0.2) g s (y) reveals a 
pronounced peak at y = la and a second, smaller one at 
y = 2cr, indicating chains of up to three particles along 
the field direction. Furthermore, we observe a ’hump’ ap¬ 
pearing directly before the first peak. This ’hump’ results 
from bended and perturbed chains as well as from a few 
clusters elongated in 2 -direction (see fig. [31(b))- Increasing 
the density the hump transforms into a pronounced peak, 
which finally overtakes the peak at y = la at p* = 0.4. 
At this density, chains in y-direction have lost their domi¬ 
nance as structural elements. Instead, clusters spread out 
in x-direction. This interpretation is supported by the 
function g s ( 2 ) which displays a similar qualitative change 
at p* = 0.4. We note that these structural details are sen¬ 
sitive to the shift parameter 8 characterizing the dipole 
positions inside the particles: Indeed, as revealed by test 
simulations, smaller values of <5 tend to surpress the lane 
formation in the ’’chaining regime” and enhance the oc¬ 
curence of chains in the ” cluster regime”. 
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Fig. 4: (Color online) Pair correlation functions between parti¬ 
cles of the same species along the x-direction (main plot) and 
the y-direction (inset). The data pertain to ff, = 80 and dif¬ 
ferent densities. In the inset, curves for different densities have 
been shifted by one a to enhance visibility. 

Mechanisms of lane formation. — In the follow¬ 
ing we concentrate on the question why the dipolar inter¬ 
actions strongly enhance lane formation compared to the 
purely repulsive SS model, as revealed by the order param¬ 
eter plots in fig. [2] We concentrate on the density p* = 0.5, 
where the differences are particularly pronounced. Be¬ 
sides the anisotropy, another main feature characterizing 
the DDP interactions is that they are effectively attractive. 
This is seen from the angle-averaged potential 

1 f 2n 

udd{t 12 ) = ^ J Udd ( r i2 , si ) S2 ) ( 2 ) 

with fi 12 being the angle of the connecting vector r 12 rel¬ 
ative to the x-axis. Numerical data for the cases si = s 2 
and Si s 2 are shown in fig. [5] Clearly both angle- 
averaged potentials display a pronounced, attractive po¬ 
tential well with nearly indentical depths. In fact, the 
angle-averaged potential somewhat resembles that of an 
LJ fluid (also plotted in fig. EJ, disregarding differences in 
the width of the attractive well and the decay of the po¬ 
tentials with ri 2 . To understand the role of the effective 
attraction in our DDP system on the lane formation we 
therefore consider, as a first step, the corresponding be¬ 
haviour of an LJ fluid. For simplicity, we choose the same 
attraction strength e* for both, different and same species. 
Results for the order parameter $ of a driven L J system as 
function of e* at fixed f£ are shown in fig. [61(a). For very 
small values of e* the order parameter is negligible, consis¬ 
tent with the behaviour of the pure SS system in fig. EJb) . 
Upon increase of e*, lanes first appear at a certain value 
and then disappear again at a substantially larger cou¬ 
pling. A similar reentrance of the non-laned (mixed) state 
occurs at other driving forces, as seen from the laning 
’’state diagram” (for p* =0.5) presented in fig. EKb). The 
state diagram moreover reveals that the onset of lane for¬ 
mation occurs at a coupling strength of e* ~ 2.5 (see dot¬ 
ted green lines in figs. EKa)-(b)) quite independent of fff 
(taking the value $ = 0.02 as a lower limit for laning), 



Fig. 5: (Color online) Angle-averaged DDP potential for si = 
S 2 and si ^ s 2 at p* = 1.58. Included is the LJ potential at 
e* = 2.3. 


while the coupling related to the disappearance of lanes 
appears to be a linear function of /)). Given the lower 
boundary e* ~ 2.5, it is interesting to make a connec¬ 
tion to the equilibrium phase diagram of the 2D LJ fluid: 
This system has a critical point (gas-liquid condensation) 
at p* = 0.35 and T* pe 0.46, corresponding to a critical 
coupling strength e* ~ 2.17 [27l[28] . The triple point (gas- 
liquid-solid) coupling strength is given as ef ps 2.5 (with 
corresponding liquid density p l f quld > 0.6 [27j). From this 
we can conclude that the non-driven LJ system at p* = 0.5 
and e* = 2.5 is in a strongly correlated state which is, in 
fact, thermodynamically unstable, i.e., it lies within the 
coexistence curve of the gas-liquid transition. Moreover, 
revisiting again fig- EKa)- (b) we see that e c provides a lower 
limit of laning. These observations suggest that the equi¬ 
librium phase separation is a prerequisite for lane forma¬ 
tion in the driven LJ system. 

We now turn to the breakdown of laning at large cou¬ 
pling strengths. Here, the external force driving two un¬ 
like particles (si sf) away from one another competes 
with the attractive LJ forces. To quantify this competi¬ 
tion we construct an effective pair interaction U^f J f(xij) 
between unlike particles, which are initially in contact 
and then driven apart in x-direction. To account for 
different vertical positions yij , we first average the x- 
component of the LJ force Flj, x = —d Xi Uhj , over the 
angle <f) between r,j and e y (at fixed distance) yield¬ 
ing F LJ:X = 7r _1 Jq sin (fFLjdf = 2 F lj /tt with F LJ = 
—d r U LJ . Integration (setting now = 0) yields the 
potential ULj{xij ) = 2/7T ULjixij). Adding the linear po¬ 
tential associated to the external force yields U^ J j:(xij) = 
ULj(xij) — fdXij. For small driving forces U^ displays a 
local minimum U m i n and maximum f7 maa; , and thus a po¬ 
tential barrier Af/ = |Umax — U m in\- Unlike particles tend 
to stick together (against the drive) if the mean kinetic 
energy per particle 2 x ksT /2 is smaller than AU(fd,e). 
Solving numerically the equation A U = ksT for various 
fl we obtain values efnaxUd) indicated by the dotted black 
line in ffg-EKb). Clearly, our simplified model describes the 
breakdown of lane formation very well. Moreover, having 
introduced , we can also understand the fact that the 
onset of laning occurs at coupling strength slightly larger 
than the critical coupling in equilibrium, e c . The driv¬ 
ing force effectively reduces the attraction between unlike 
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Fig. 6: (Color online) (a) Laning parameter versus coupling strength for the LJ system at // = 7, p* = 0.5. Inset: Enlarged 
view for small e*. (b) Laning ’’state diagram” of the LJ system. In (a) and (b), red, dotted green and black lines correspond 
to e*, tl, and respectively, (c) Corresponding state diagram of the DDP system. The dotted white, red(blac.k) lines 

correspond to p* and at yo = 0.645(0.9)cr . 


particles, yielding a shift of phase separation. 

Given this background, we now consider in fig. 0c) the 
state diagram of the full, anisotropic DDP fluid (p* = 
0.5) in the parameter plane spanned by driving force and 
coupling strength. The latter is now given by p* (the 
soft-sphere repulsion is fixed). Clearly, the structure of 
the state diagram resembles that of the LJ system. In 
particular, the onset of laning depends only weakly on 
and there is an upper boundary for p* beyond which 
laning disappears. 

This prompts the question whether the onset of lan¬ 
ing in the driven DDP system can be related to an equi¬ 
librium fluid-fluid phase separation, as in the LJ sys¬ 
tem. We note that already the equilibrium DDP sys¬ 
tem is a true binary mixture, since the full pair inter¬ 
actions depend on s», see fig. 0 Therefore the possi¬ 
ble fluid-fluid transitions (if existent at all) are, in gen¬ 
eral, combinations of condensation and demixing. In the 
present study we did not carry out simulations to ex¬ 
plore these questions properly. Still, one can perform 
some estimates based on mean-field density functional the¬ 
ory (DFT). We focus on the occurence of condensation 
(rather than demixing), because the angle-averaged (i.e. 
mean-field) interaction between unlike DDPs is nearly the 
same as that between like particles (see fig. 0. The cor¬ 
responding stability condition of the homogeneous, mixed 
phase is that the isothermal compressibility, xt, has to 
be positive. According to Kirkwood-Buff theory [29] 
one has (for a symmetrix binary mixture composed of 
species A and B) x^ 1 oc 1 - (p/2 ) (cUa(O) + cUb(O)), 
where caa(ab) (0) are the Fourier transforms of the di¬ 
rect correlation functions (DCFs) CAA{AB){ r 12) in the 
limit of long-wavelengths (k -A- 0). In our case, AA 
( AB ) corresponds to si = S2 (si = —s 2 ). Further¬ 
more, we approximate the DCFs according to a random 
phase approximation, that is, c SlS2 (r 12) = CHs{ r 12 )^(<7 — 
ri 2 )-{k B T)~ 1 U DD (r 12l si, s 2 )0(r 12 -a), where c H s is the 
Percus-Yevick DCF of a pure hard-sphere fluid [3D]. The 
Fourier transforms of the second, mean-field like contribu¬ 
tion to the DCFs yield essentially the spatial integral over 


the effective potentials defined in eq. 0- Numerical inves¬ 
tigation of the resulting expression for xt at p* = 0.5 re¬ 
veals that, upon increasing p* from zero (where \t reduces 
to the hard-sphere compressibility), xt becomes indeed 
negative at p* = 1.31, indicating an instability (’’spinodal 
point”) related to condensation. Thus, simple mean-field 
DFT predicts the existence of a condensation transition in 
the equilibrium DDP fluid. Further, for p* = 0.5, the tran¬ 
sition should occur at a coexistence value p* oex < p*\ its 
actual value is, however, more eloborate to determine [3Tj . 
BD test simulations of the non-driven system at p* = 0.5 
and dipole moments larger than p* reveal indeed a phase 
separated structure consisting of thick columns (involving 
both particle species) and large voids in between. 

The question now is, does p* play a decisive role for 
the driven DDP system? Considering fig. 0c) we find 
that p* (indicated by a white line) yields indeed a good 
estimate for the onset of laning at low values of the driving 
force < 30). This suggests that, similar to the LJ 
system, equilibrium phase separation is a prerequisite of 
laning at the density considered. Only for larger driving 
forces (/^ > 30) the estimate worsens (this effect was not 
observed in the LJ system where, however, calculations 
where restricted to small /£). 

The breakdown of laning at large p* can be esti¬ 
mated, similar to the LJ system, by constructing an one- 
dimensional (x-dependent) potential, in which the DDP 
interaction between unlike particles competes with the 
driving force. We consider examplary configurations yij = 
yo (a simple force average over different configurations yij 
is not appropriate due to the interaction anisotropy) and 
obtain the potential U^ p (xij) = U B Dp(x',yo) — fd^ij 
which displays a barrier AU(yo 1 fd,p)- Figure 0c) in¬ 
cludes the numerical solution Pmax(fd) (dotted red line) 
of the equation A U = k B T at yo = 0.645cr, correspond¬ 
ing to the most attractive, ’’staggered” configuration at 
contact (see Fig. 0b)). Clearly, this line describes the 
breakdown qualitatively, but overestimates the influence 
of attraction. Quantitatively better results are obtained 
with yo = 0.9ct (dotted black line in fig. 0c)), correspond- 
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ing to a less attractive initial configuration. The ’’true” 
potential should be seen as a (weighted) average over dif¬ 
ferent yo. Finally, the reduction of attraction expressed 
in the effective potential explains the shift of the onset of 
laning with larger 

Conclusions. — Our study reveals new und unex¬ 
pected aspects of laning in driven colloidal mixtures with 
anisotropic, dipolar interactions between the particles. 
Most importantly, anisotropic interactions do not sur- 
press the lane formation despite the strong tendency of 
the particles to self-assemble into chains (which, in the 
present model, are oriented perpendicular to the force). 
Rather, we observe stabilization of laning (relative to re¬ 
pulsive systems) in a window of interaction strength, with 
the anisotropic interactions influencing strongly the struc¬ 
ture of the lanes. This stabilization is due to the ef¬ 
fective, angle-averaged attraction inherent in our model, 
which eventually leads to equilibrium phase separation 
(condensation). Interestingly, the corresponding coupling 
strengths also give a good estimate for the onset of lan¬ 
ing upon at fixed driving force, and we have observed 
similar behaviour in the simpler (isotropic) LJ system. 
This indicates that laning in systems with attractive in¬ 
teractions is intimately related to phase separation, irre¬ 
spective of their detailed form. Our observations should 
be directly measurable in real systems of field-propelled 
metallodielectric particles [l5j. Moreover, we expect our 
findings to be transferable to other driven colloidal sys¬ 
tems with direction-dependent interactions, examples be¬ 
ing ’’patchy” particles which display both equilibrium ag¬ 
gregation and condensation |32l , or dipolar colloids under 
shear flow [331 . 

Starting from the present study, there is several fur¬ 
ther intriguing questions to be explored. For example, 
is the behaviour observed here related to a Raleigh-Taylor 
(RT) instability occuring in driven, macroscopically phase- 
separated mixtures? In the context of RT, the criti¬ 
cal wavelength [34] can be estimated by a simple crite¬ 
rion involving the surface tension; this might provide a 
way to quantify the lane width of our system.Indeed, the 
present simulation data are not yet conclusive to predict, 
e.g., a relation between lane width and dipolar coupling. 
We also point out possible connections to driven,phase- 
separating magnetic mixtures [J5]. Finally it would be 
worth to explore connections to active fluids consisting of 
self-propelled agents with attractive interactions. Indeed, 
the interplay of self-organization in such systems and equi¬ 
librium phase separation is currently a very lively field of 
research [T9lf20lf36Ff38] . 
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